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Abstract 

We address the consistency of Hofava's proposal for a theory of quantum gravity 
from the low-energy perspective. We uncover the additional scalar degree of freedom 
arising from the explicit breaking of the general covariance and study its properties. 
The analysis is performed both in the original formulation of the theory and in the 
Stiickelberg picture. A peculiarity of the new mode is that it satisfies an equation of 
motion that is of first order in time derivatives. At linear level the mode is manifest only 
around spatially inhomogeneous and time-dependent backgrounds. We find two serious 
problems associated with this mode. First, the mode develops very fast exponential 
instabilities at short distances. Second, it becomes strongly coupled at an extremely low 
cutoff scale. We also discuss the "projectable" version of Hoi^ava's proposal and argue 
that this version can be understood as a certain limit of the ghost condensate model. 
The theory is still problematic since the additional field generically forms caustics and, 
again, has a very low strong coupling scale. We clarify some subtleties that arise in the 
application of the Stiickelberg formalism to Hofava's model due to its non-relativistic 
nature. 
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1 Introduction and summary 



Recently, Horava has proposed a new approach to the theory of quantum gravity [T]. The 
key idea of the proposal is to equip space-time with a new structure: a foliation by space- 
like surfaces. This foliation defines the splitting of the coordinates into "space" and "time" 
and breaks the general covariance of general relativity (GR). Then one can improve the 
UV behavior of the graviton propagator and ultimately make the theory power-counting 
renormalizable by adding to the GR action terms with higher spatial derivatives. At the 
same time the action in the ADM formalism contains only first order time derivatives, which 
allows to circumvent the problems with the ghosts appearing in covariant higher order gravity 
theories [2] . The higher derivative terms naively become irrelevant in the infrared and it was 
argued in [1] that the theory reduces to GR at large distances. 

However, the consistency of the above proposal is far from being clear. The main concern 
comes from the fact that the introduction of a preferred foliation explicitly breaks the gauge 
group of GR down to the group of space-time diffeomorphisms preserving this foliation. As 
already pointed out in [1] this breaking is expected to introduce extra degrees of freedom 
compared to GR. The new degrees of freedom can persist down to the infrared and lead 
to various pathologies (instabilities, strong coupling) that may invalidate the theory. An 
illustration of this phenomenon is provided by theories of massive gravity where special care 
is needed to make the additional degrees of freedom well-behaved [31 HI [5]. 

In the recent works in the topic there have been several controversial claims about the 
properties of the extra degrees of freedom. In [6] the new mode was identified among the 
perturbations around a static spatially homogeneous background in the presence of matter. 
The mode was argued to be strongly coupled to matter in the limit when the theory is 
expected to approach GR, making it hard to believe that a GR limit exists. It is worth noting 
that the mode found in [6] is not propagating: its equation of motion does not contain time 
derivatives [7]. Thus it remains unclear from this analysis whether this mode corresponds 
to a real degree of freedom or can be integrated out as unphysical. The observation that the 
extra mode is non-propagating was generalized in [8] to the case of cosmological backgrounds. 
The interpretation of this result given in [8] is that actually the Hofava gravity is free from 
additional degrees of freedom. It was also claimed that the strong coupling is alleviated by 
the expansion of the Universe. Finally, the non-linear Hamiltonian analysis performed in 
[9] shows that the phase space of Hofava gravity is 5-dimensional. This result is puzzling: 
a normal degree of freedom corresponds to a 2-dimensional phase space; so the result of [9] 
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suggests that the number of degrees of freedom in Hofava gravity is two and a half. Two of 
these degrees of freedom are naturally identified with the two helicities of graviton. But the 
physical meaning of the extra "half-mode" is obscure. 

The aim of the present paper is to clarify this issue. We show that Hofava gravity does 
possess an additional light scalar mode. For a general background the equation of motion of 
this mode contains time derivatives implying that the mode is propagating. The peculiarity 
of Horava gravity is that the equation for the extra mode is first order in time derivatives. 
Still, the solution corresponds to waves with a background dependent dispersion relation and 
is fixed once a single function of spatial coordinates is determined as the initial condition in 
the Cauchy problem. This explains why this mode corresponds to a single direction in the 
phase space. 

Next we address the consistency of the Hofava proposal by studying the infrared proper- 
ties of the extra mode. We find that its dynamics exhibits a number of bad features. First, 
the mode becomes singular for static or spatially homogeneous backgrounds. Namely, the 
mode frequency diverges in that limit. This explains why this mode has been overlooked 
in the previous analyses of perturbations in Hofava gravity P, El IB]- Second, for certain 
(background-dependent) values of spatial momentum the mode becomes unstable. Again, 
the rate of the instability diverges if one takes the static / spatially homogeneous limit for 
the background metric. Third, we show that at energies above a certain scale the extra 
mode is strongly coupled to itself, and not only to matter. We find that the strong coupling 
scale is background dependent and goes to zero for flat / cosmological backgrounds. Hence, 
the model suffers from a much more severe strong coupling problem than pointed out in [6j, 
where the dependence of the strong coupling scale on the background curvature was ignored. 
Because of the strong coupling the Hofava model can be trusted only in a narrow window 
of very small energies, way below the Planck scale; this point is illustrated schematically in 
Fig. [H This implies that the Hofava model cannot be considered as consistent theory of 
quantum gravity. 

To unveil the properties of the extra mode we make use of the Stiickelberg formalism. For 
the case at hand the Stiickelberg trick is synonymous to the covariantization of the model. 
As a result we obtain a scalar-tensor theory with the time derivative of the scalar field 
developing non-zero vacuum expectation value. The invariance under foliation preserving 
diffeomorphisms implies that the theory has an internal symmetry consisting in reparameter- 
izations of the scalar. We clarify the subtleties that arise in the application of the Stiickelberg 
procedure to Hofava gravity due to the intrinsically non-relativistic nature of the proposal. 
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Figure 1: The energy scales appearing in the Hofava modeh the Planck mass Mp and the 
strong coupling scale Asc- The theory can be trusted only in the narrow window of energies 
< E < Asc (dashed region). The scale A^c is much smaller than Mp, it depends on the 
background curvature and goes to zero for flat / cosmological backgrounds, see Eq. fl5^ 
below. 

The covariantization of the higher space derivatives of the model leads to higher covariant 
derivative operators in the equations of motion. Naively, this would imply the appearance of 
too many degrees of freedom. However, in the theory at hand the higher derivative opera- 
tors are of a special type that allows for a well-posed Cauchy problem with reduced number 
of initial data in the preferred foliation. In this way the number of degrees of freedom is 
decreased and matches with the number of modes in the non-covariant formulation. As a 
byproduct, we point out a large class of covariant higher derivative operators that allow for 
the reduction of the number of degrees of freedom in a preferred Cauchy slicing. 

There exist two versions of the original Horava proposal [Tj. The difference between 
them lies in an additional restriction which can be imposed on the lapse function. Namely, 
one can require the lapse to be "projectable" , i.e. be constant along the foliation surfaces. 
In the present paper we are mainly interested in the non-project able case; all the previous 
discussion refers to this case. The projectable version of the theory is briefly discussed at the 
end of the paper. We argue that in this case the Hofava gravity is equivalent to a specific 
limit of the ghost condensate model [10]. This implies that now the theory possesses a full- 
fledged extra scalar degree of freedom with second order equation of motion. The classical 
dynamics of linear perturbations of the scalar is regular for all backgrounds. In this sense the 
projectable version of the theory is better behaved than the non-projectable one. However, 
beyond the linear or classical level the additional scalar still exhibits pathologies. As shown 
in [11], at the classical level the dynamics of the projectable version of Hofava gravity is 
equivalent to GR supplemented by a pressureless fluid. As we discuss, the fluid component 
is precisely described by the extra scalar, the fluid velocity being proportional to the scalar 
gradient. A well-known property of pressureless fluid is to develop caustics where the fluid 
velocity becomes ill-defined. For the scalar at hand it means that the theory inevitably 
breaks down after finite amount of time evolution. At the quantum level, the extra mode 
exhibits unacceptably low scale of strong coupling. We comment on the possible ways to 
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address these problems. 

The paper is organized as follows. In Sec. [2] we describe the model and formulate the 
Cauchy problem for it. In Sec. [3] we derive the linearized equations for perturbations about 
an arbitrary background. We find an explicit expression for the extra mode and show that 
it obeys first order equation in time. In Sec. H] we turn to the Stiickelberg analysis of the 
model which allows us to study the properties of the extra mode in a transparent way. We 
discuss the subtleties in the application of the Stiickelberg formalism to the case at hand. 
In Sec. |5] we take the limit when the GR part decouples from the Stiickelberg sector and 
concentrate on the latter. We find that the linear perturbations of the field exhibit fast 
exponential instability and unacceptably low scale of strong coupling. Finally, in Sec. M we 
consider the projectable version of Hofava gravity and argue that it is equivalent to a specific 
limit of ghost condensate model. Sec. [7| contains concluding remarks and discussion of future 
directions. Technical details are deferred to Appendix lAl 



2 Cauchy problem for Hofava gravity 

We consider the class of non-relativistic generalizations of GR proposed in [T]. One starts 
with the ADM decomposition of the space-time metric, 

ds^ = (iV^ - NiN')dt^ - 2Nidx'dt - -fijdx'dx^ . 

Then the action for a theory of this class can be written in the generic form 
/Vf2 r 

S=^l d'^v^ AT {Ki,K'^ - + iR + CR,,R'^ + ...), (1) 

where Mp is the Planck mass; Kij is the extrinsic curvature tensor for the surfaces of constant 
time, K is its trace; Rij, R, 7 are the Ricci tensor, Ricci scalar and the determinant of the 
spatial metric 7^^, and A^ is the lapse function. The extrinsic curvature is related to the time 
derivative of the metric in the usual way, 

= ^ - V^N, - V,N,) . (2) 

Throughout the paper, if not stated otherwise, 3-dimensional indices . . . are raised and 
lowered using 'jij, and the covariant derivatives carrying these indices are understood in the 
3-dimensional sense. The ellipsis in ([T]) represents higher order terms constructed out of the 
metric 'jij using only spatial derivatives and invariant under 3-dimensional diffeomorphisms. 
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As discussed in [Ij, the introduction of terms of sufficiently high order (with six spatial 
derivatives) yields a theory which is naively power-counting renormalizable. 

The purpose of the present paper is to study the properties of the theory ([!]) in the 
infrared. The precise structure of the higher order terms is not essential for this analysis. To 
illustrate the general effect of these terms it suffices to consider explicitly one of them which 
we choose to be the square of the Ricci tensor. 

The action ([T]) reduces to GR for the values of the parameterj^ \ = I, ( = 0. Away 
from these values ([T]) explicitly breaks general covariance down to the subgroup consisting 
of spatial diffeomorphisms and reparameterizations of time, 

x(t,x) , t i{t) . (3) 

These transformations preserve the foliation of the space-time by the surfaces t = const. Note 
that if one assumes that the action ([T]) defines a consistent quantum theory, one expects all 
the parameters in this action to acquire radiative corrections and flow by the renormalization 
group. In particular, the parameter A is expected to be generally different from 1. It was 
argued [I] that in the infrared limit the higher order terms become negligible and so one 
may expect to recover GR if the parameter A flows to 1 in that limit. We will see below 
that this expectation is incorrect: the explicit breaking of general covariance leads to the 
appearance of an extra degree of freedom in the infrared which becomes strongly couplecj^ 
when A approaches 1. 

Varying the action with respect to A^, Ni, ■yij yields the following equations, 

- KijK'^ + XK^ + ^R + CRijR'' = , (4) 

ViK'^ - XV^K = , (5) 
d 

- - Xj'^K) - (1 - 2X)NKK'^ - 2NK'^Kl + NK^K'^^Y^ - ^NR'^ 

- ^Y^AN + ^V'V^N - 2CNR''Ri - C"^ i{N R''^)Y^ - CA{NR^) 

+ C[VkV'{NR''^) + VkV^iNR''')] = . (6) 

Here A = 7*-' VjV,- and we have flxed the gauge Ni = 0. This system has to be supplemented 



^The value of the parameter ^ is not important: it can always be absorbed into the rescaling of the time 

coordinate and the shift Ni. 

^More precisely, the additional mode is weakly coupled only in a narrow window at low energies. This 
window depends both on A and the parameters of the background geometry; it shrinks to zero both when 
A ^ 1 or when the background curvature vanishes, see Sec. El 
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by the evolution equation ([2]) for the metric which in the chosen gauge takes the form 

= 2NKij . (7) 

Let us analyze the Cauchy problem for the system (jlj) - (jTj). The set of initial data at 
t = consists of the values for 7jj, K^^ and at this time. The initial data have to 
satisfy the constraints (jl]) and ([5]). Then Eqs. ([6]), ([Tj) describe the time evolution of the 
extrinsic curvature and the metric. However, the system (jlj) - ([7]) is incomplete: it does 
not allow to determine the time evolution of the lapse A^. In GR this ambiguity is a gauge 
artifact removed by appropriate gauge fixing, e.g. = 1. In our case the gauge freedom 
is absent, and is a genuine dynamical field. To obtain the missing equation one notices 
that, due to the lack of gauge invariance, the Hamiltonian constraint (HI) is not automatically 
preserved by the time evolution. Imposing that the constraint holds at any time one gets a 
secondary constraint which produces an equation for A^. Taking the time derivative of (jl]) 
and simplifying the result with the use of the rest of the equations we obtain 

- KkjV'^R'^ + R'^^VkK} - K'NkR^' + R^V^Kkj)] 1=0. 



As expected, the l.h.s. of this equation vanishes identically in the case of GR as a consequence 
of gauge invariance. For A and C away from their GR values this equation allows to determine 
the lapse at any moment of time provided the configuration of 7jj, K"^^ is given; in this way it 
imposes additional constraint on the initial data. It is important to notice that the constraint 
vanishes whenever the extrinsic curvature or gradients of the curvature tensors are zero. In 
particular, this happens for spatially homogeneous or static configurations. Note also that 
Eq. ([8]) has the form of the conservation of a (space-like) current. This fact acquires a 
natural interpretation in the covariant picture where (IE]) becomes the equation of motion 
of the Stiickelberg field, and the current corresponds to the shift symmetry of this field, cf. 
Eq. (Hg) below. 

The system (H]) - ([8]) constitutes the complete set of equations of motion for Horava 
gravity. Let us count the number of independent Cauchy data for this system. Originally 
the set of initial data for 7jj, K^^ ^ N contains 6 + 6 + 1 = 13 functions of spatial variables x. 
The constraints (|1D, ([5]), ([8]) eliminate 1 + 3 + 1 = 5 of them. Additionally, 3 functions are 
removed by the residual (time-independent) gauge transformations of spatial coordinates. 
Thus we are left with 13 — 5 — 3 = 5 arbitrary functions as initial data. 4 of these functions 
are identified as initial data for the two helicities of the graviton. The remaining freedom 
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in the choice of one more function imphes the presence of an extra mode which is absent 
in GR. Note that because the Cauchy data for this mode are hmited to a single function, 
the corresponding evolution equation must be first order in time. This agrees with the 
observation made in [9J that the phase space of the Hofava gravity has odd dimensionality. 
Our task below is to investigate the properties of the extra mode. 

3 The elusive mode 

In this section we reveal the extra mode explicitly. For the sake of the argument we restrict 
to the case when the higher order terms in the action are absent, ( = 0, and breaking of 
general covariance arises only from A 7^ 1. This restriction allows to capture the essential 
physics of the extra mode, without overloading the paper with lengthy formulae. As it will 
become transparent below, the main conclusions are unaffected by the values of ( and A, so 
long as the general covariance is broken by some of these terms. For C = the secondary 
constraint ([H]) reduces to 

Vi {N'^V'K) = 0. (9) 

Consider small perturbations of the fields about a background ■yij,K^^,N. We assume that 
the background satisfies the equations of motion of the Hofava gravity, but is arbitrary 
otherwise. Thus we write 

N = N + n . 

The next step is to plug these expressions into Eqs. dl]) - ([8]) and expand them to linear 
order in the perturbations hij, k,^\ n. We obtain 

- 2KijK'^ - 2KiK^^hij + 2\Kk + 2XKK'^hij - ^Ah + ^V'V^hij - ^R'^hij = , (10) 
ViK^' -XV'K + K''^Wkhi-^-^K''^V^hki + ^K''^Wkh + XV^Khi = , (11) 
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(1 - 2X)NKk'^ - 2NKiK^'j - 2NKiK^' - iV((l - 2\)K'^ - 2Xf^K)K 
XNK^h'^ - 2NK'^K^^hki - iV((l - 2X)K'^ - 2Xf^K)K^^hki 



+ AiV + V'^ V'iV hki - V^V-^'iV h\ - V^V'iV 



, (12) 



hij = 2NKij + 2NK^hik + 2NK^hjk + 2Kijn , (13) 
2ViKV'n + AKn + 2V,iVVV + N Ak + NK'^Ahij + 2NV^K'^Vkh^j - NVjKVih'^ 

- 2V'NV^Khij + 2V'NViK^^hki = . (14) 

Here the indices are raised and lowered using the background metric 'jij, and the covariant 
derivatives are understood with respect to this metric. 

One makes an important observation. As discussed above, Eq. ( |T4|) is supposed to deter- 
mine the evolution of the lapse. However, the terms linear in n disappear from this equation 
whenever the gradients of the background extrinsic curvature vanish. In particular, this 
happens for static or spatially homogeneous backgrounds. Then, instead of determining the 
lapse, Eq. (IT^ imposes a constraint on an otherwise propagating field, making the extra 
mode non-dynamicaH 

It is instructive to work out the case of Minkowski background in a certain detail. In this 
case Eq. ( IT2|) takes the form 
d_ 
di' 



didjU — 5ijAn 

(15) 

. 



+ ^Ahij - ^dkdihjk - ^dkdjhik + ^didjh - SijAh + Sijdkdihki 



The fact that the extra mode does not propagate at hnear level in homogeneous or static backgrounds 
also holds for the general Hof ava action ([1]) . The reason is that the combination in the square brackets in the 
secondary constraint ([5]) vanishes on these backgrounds. This is precisely the combination which multiplies 
the perturbation of the lapse in the linearized equation. 
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The linearized Hamiltonian constraint (ITOl) yields A/i = didjhij; therefore the trace of (ITSl) 
reads 

(1 - 3X)k = -2^An. 

If n were determined by dHl), this would be a first order equation for k. However, in 
Minkoswki (HM reduces to 

Ak = , 

which restricts the perturbation of both the extrinsic curvature and the lapse to vanish. The 
rest of the argument proceeds in the same way as in GR and one concludes that all scalar 
modes are non-propagatinj^. The same effect occurs in any spatially homogeneous or static 
background. This explains why the extra mode was overlooked in the previous analyses 
[SI [71 E] that focused on this class of backgrounds. 

According to the above discussion the extra mode reveals itself only in backgrounds which 
are both time-dependent and spatially inhomogeneous. Finding an exact solution of Hof ava 
gravity with these properties is a difficult task. Fortunately, for our purposes it is not needed: 
it is enough to realize that such backgrounds exist. As a concrete example, one can keep in 
mind a large gravitational wav^. 

In the generic case, the system of linearized equations is intractable. What saves the day 
is the fact that the extra mode appears in any background such that the terms proportional 
to n in ( fT4|) do not cancel. It is enough to consider perturbations at space-time scales much 
shorter than the characteristic distance of the variation of the background. This allows to 
treat the background fields as almost constant at the scales of interest and ensures the validity 
of Fourier analysis at these scales. Technically this amounts to keeping in the equations only 
terms with least number of derivatives of the background. Let us make this point more 
quantitative. We assume that the background metric changes at characteristic space-time 
scale L. Then we have Rij ~ ^/L"^, Kij ~ 1/L. We are interested in perturbations at 
distances much shorter than L. This means that we consider perturbations with frequencies 
and momenta uj,p ^ 1/L. Consequently, in Eqs. ( JTOl) - (|T4l) we can neglect terms with 
derivatives of the background in comparison to the terms with derivatives of perturbations. 
This does not imply throwing away all the terms with background gradients: some of these 

"^This statement is true only at linear order in perturbations. The non-linear corrections will bring back 
the propagating mode as it is clear from the study of perturbations in general backgrounds and from the 
Stiickelberg picture, see below. 

^Moreover, once the perturbations to any metric are considered, the presence of inhomogeneities is uni- 
versal. 
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terms may be the leading ones. For instance, the first term in Eq. (fT4|) is the leading 
contribution containing the lapse n in this equation. In this way we obtain the simplified 

system, 

- 2KijK'^ + 2XKk - ^Ah + id'd^hij = , (16) 
diK'^ - + K^^'dkh^i - ^-^K^'d^hki + ^K'^dkh = , (17) 



+ —Ah'^ - —d'^d'hi - —d'^d^hl + —d'd^h - NS'^Ah + N5'^d^d^hki 
2 2 2 2 

hij = 2NKij + 2Kijn , (19) 

2d'nViK + NAk + NK'^Ahij = , (20) 



d'd^n - 5'^ An 



where without loss of generality we have seto Sij. Let us first consider Eq. f|T9l) . It 

follows from this equation that either or Kij-n is at least of order uhij. Let us further 
assumjl uj ^ p: we will see shortly that the dispersion relation of the extra mode obeys 
this inequality. Then we see that in Eq. (l20l) the last term is always negligible. A similar 
reasoning shows that one can neglect all the terms containing hij in (fT8|) . This yields a closed 
system of equations for k*-' and n: 

-k'^ + XS'^k + e [d'd^n - S'^An] = , (21) 
2d'nViK + NAk = . (22) 

We point out that this system is explicitly first order in time derivatives. As already men- 
tioned, at short scales we can treat the combinations of the background fields appearing in 
the above equations as constant. One performs the Fourier decomposition 

/iij-,K^^noce-*'^*+*P^ 

and finds that the solution of the system (^T^ - fl22]) has frequency 

= 1 ^ — - . 23 

(1 - 3X)p^djK ^ ' 



^The background metric can be always brought to this form in the vicinity of any given point by the 
time-independent 3-dimensional diffeomorphism. 

^The analysis in the complementary regime w ~ p reveals only the two transverse traceless modes of the 
graviton. 
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The extrinsic curvature for this solution is determined in terms of the lapse, 

'^'^ = '-^^ ( - (1 - 3A)py + (1 - A)5V)n . (24) 

Note that uo ~ p^L!^ which is indeed much bigger than p. Besides, the extrinsic curvature 
behaves as 

~ . (25) 

From this estimate one concludes that the r.h.s. of 0191) is dominated by the second term. 
This yields for the perturbations of the metric 

2i(l - 2>\)p^dkK 

Finally, we have to check the constraints flTBl) . f[T7|) . From fl2^ we see that 

n 



hij = K,,n . (26) 



hi 



{pLf 

This estimate together with (l25l) implies that the terms containing hij in (fT7j) should be 
neglected compared to the first two terms; on the other hand, all the term in f|T6|) are of the 
same order. Then it is straightforward to verify using the explicit expressions ^2M . ( 126|) that 
the constraints are satisfied. 

Let us briefly summarize our results. Eqs. flM|) and fl26|) provide the explicit expression for 
the extra mode of Hofava gravity in the short wavelength limit. This mode is parametrized 
by a single scalar function n(p) and has the dispersion relation (!23|) . Together with the two 
polarizations of the graviton found in the complementary regime u; ~ p this matches with 
our counting of degrees of freedom in Sec. O The frequency ( l23i) of the mode diverges when 
the gradients of the background extrinsic curvature vanish. Thus the mode becomes singular 
for spatially homogeneous or static backgrounds. 

The expression (1231) also diverges for the modes with spatial momenta perpendicular to 
the gradient diK. Naively, one could try to find the behaviour of these modes by applying 
the Fourier analysis to the system (fT6l) -( l20ll including next to leading term from Eq. ([HI), i.e. 
the term AKn. However, this would give an incorrect result. The leading-order expression 
(!23|) for the frequency of the mode depends on the background fields and hence on the space- 
time point where one performs the Fourier decomposition. Therefore at the subleading level 
the first term in Eq. ( fT4l) produces contributions of the form 

ntdiCud'^K ~ nutL'^ , 
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which at the time scales of interest t ~ are of the same order as the term 



AKn ~ nL"^ . 

Thus, the consistent treatment of the subleading effects requires going beyond Fourier anal- 
ysis. Instead, one has to implement the WKB expansion in order to properly account for 
the inhomogeneity of the background. We will perform this study in the decoupling limit in 
Sec. El where we will find that the subleading corrections generically lead to fast exponential 
instability of the extra mode. 

To close this section, let us emphasize that the qualitative properties of the extra mode 
mentioned above are generic for any Hofava-type Lagrangian, even though in our consider- 
ation we mainly focused on the term (1 — X)K^. Technically, the important point is that 
the perturbation of the lapse n enters in the linearization of the secondary constraint (8) 
with a coefficient that does not vanish on inhomogeneous and time- dependent backgrounds. 
Then this equation relates n with the perturbation of the trace of the extrinsic curvature k 
and leads to a first order evolution equation. In other words there will be always analogs 
of Eqs. i^I^, fl22p whenever the Lagrangian contains terms violating the four- dimensional 
general covariance. This is what makes the scalar mode propagate at linear level around 
inhomogeneous and time-dependent backgrounds for any Horava-type Lagrangian. 



4 Stiickelberg formalism 

To get more insight into the dynamics of Horava gravity we use the Stiickelberg formalism. 
This will allow us to clearly separate the extra propagating mode and perform a detailed 
analysis of its properties. 

The first step is to restore the full general covariancqj at the expense of introducing the 
corresponding Stiickelberg field. Namely, we encode the foliation structure of Hofava gravity 
in a scalar field 0(a;) with non-vanishing time-like gradient. The surfaces of the foliation are 
then defined by the equations 

= const . (27) 

The original action ([1]) is written in the gauge where the field coincides with time, (p = t. 
Below we will refer to this choice of coordinates as "unitary gauge" . 



^While this paper was in preparation Ref. [12] appeared which also deals with the topic of the covariant 
form of the Hofava gravity. 
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Before obtaining the explicit expression for the action in an arbitrary gauge let us antici- 
pate some of its properties. First, due to the presence of the new field 0, we expect the action 
to be some kind of tensor-scalar theory. Second, the invariance of the original formulation 
([T]) under time reparameterizations (second equation in Q) translates into the symmetry of 
the covariant action with respect to reparameterizations of the Stiickelberg field, 

= , (28) 

where / is an arbitrary monotonous function. The appearance of a time- dependent vev for 
breaks the product of this symmetry and general covariance down to the diagonal subgroup. 
The latter translates in the unitary gauge into the invariance under foliation preserving 
diffeomorphisms ([3]). 

To proceed one notices that the quantities appearing in ([T]) are the standard geometrical 
objects (induced metric, extrinsic and intrinsic curvature) characterizing the embedding of 
the hypersurfaces defined by fl27|) in space-time. The central object in the construction of 
these quantities is the unit normal vectoj^ u^. Explicitly, 

Un = —7= 5 

where 

Note that is automatically invariant under the transformations fl25]) . Other geometrical 
quantities associated to the foliation are constructed out of and its derivatives. We have 
the following expressions for the spatial projector: 

the extrinsic curvature: 

/C = P Vu = ^ pf'P^XI V 6 



and the intrinsic Riemann tensor: 

7^ V = n P'p p! ^'^R'p.s + ^^u. - /C^^/C.p , (29) 



^Throughout the paper the Greek indices . . . are raised and lowered using the 4-dimensional metric 
Qpu while the Latin indices i, j, . . . are raised and lowerd using the spatial metric jij . The same correspondence 
applies to the covariant derivatives carrying these indices. 
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where in the last equation ^'^'>R°'p^^ is the 4-dimensional Riemann tensor. Now it is straight- 
forward to obtain the covariant form of the action ([T]) by identifying the quantities appearing 
in the ADM decomposition with the appropriate combinations of u^, P^y, KL^u, etc. in the 
unitary gauge. For instance, in this gauge one has 

Mo = ^ = iV , Mi = , (30) 
_ n p*-7 — —o/*-? 



pOO _ pOi _ g pij _ 



In this way from ([T]) we obtain the following covariant action, 

s = ^J d'x^[ - (^)i? + (1 - A)/c2 + cp^'^pf-n^.n,^ + ...} 



(31) 



where )C = JC^, and 



n,, = p^p^ (^)i?,^ - u^uf" + /c^/c^p - /c^/c,p . 

In deriving fl3Tl) we set for simplicity ^ = 1; as we have already mentioned the value of this 
parameter is not physically relevant and we will stick to this choice from now on. The above 
action describes gravity interacting with a derivatively coupled scalaiF^ 0, which enters into 
( 13T|) through the combinations P^^^, IC^^, Tl^v 

From the action ( 13T|) the advantage of the Stiickelberg formalism is clear: it allows to 
transfer the extra mode of Horava gravity from the metric sector to the 0-sector. Indeed, 
due to the general covariance of the action one can always choose the gauge where the metric 
sector contains only the two transverse traceless modes of the graviton. At the same time 
the extra mode is unambiguously identified with the fluctuation of the foliation structure. 

At this point we encounter a puzzle. To be consistent with the counting of degrees of 
freedom in the unitary gauge the equation of motion for the Stiickelberg field must be first 
order in time derivatives. On the other hand, it is easy to see that the action (13T|) contains 
more than two derivatives of the field (f). For example, consider the term proportional to 
(1 — A). Written explicitly in terms of it reads. 



Ml(l-X) f,,^^l V^0V-0 



= I d^xV^ ^ ( 00 - ^^^V,V.0 ) , (32) 



Clearly, the terms in the action containing the scalar disappear when A = 1 and all the higher order 
terms vanish, i.e. in the pm'e GR case. 
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and contains four derivatives! I Thus, for general choice of space-time coordinates, the equa- 
tion of motion for is fourth order in time derivatives. The corresponding Cauchy problem 
requires four arbitrary initial data; this apparently contradicts the counting of degrees of 
freedom performed in Sec. [21 where we found only one additional function compared to GR. 
The resolution of this puzzle lies in the fact that the higher- derivative equation following 
from fl32|) is of a very special type. There exist a particular choice of coordinates for the 
formulation of the Cauchy problem where less initial data are required. This Cauchy slicing 
is precisely the preferred foliation of the model. In these coordinates the number of time 
derivatives in the equation for (p is reduced to one, which matches with the counting of the 
degrees of freedom in the unitary gauge. 

To illustrate the point about the reduction of degrees of freedom in a preferred frame let 
us make a digression and consider the following equation for a non-relativistic scalar ip, 

<^ + (_l)9AV = 0, (33) 

where g > 2 is an integer number. This equations describes one degree of freedom and 
the corresponding Cauchy problem involves two initial data. The general solution of the 
equation (l33l) is a collection of waves with dispersion relation tu^ = p^'^. 

However, when one attempts to write down the same theory in a manifestly Lorentz 
invariant way one seemingly encounters the problem that the theory is higher derivative 
both in space and time directions. Indeed, in a generic Lorentz frame Eq. (1331) acquires up 
to 2q time derivatives of (p. An observer in this frame would conclude that the number of 
degrees of freedom is g > 1, as to solve the equations of motion she would need up to 2q 
initial conditions. Does this contradict the counting of degrees of freedom in the original 
frame? The answer is no, because the two formulations of the Cauchy problem are physically 
inequivalent. Indeed, recall that specification of the solutions which are considered as physical 
involve fixing the boundary conditions at spatial infinity. In the simple case of Eq. fl33l) the 
natural choice is to impose vanishing of the field, 

— > , at |x| ^ oo . (34) 

When one thinks of the Cauchy problem in the boosted frame, one also implicitly assumes 
the condition (IM]) but now imposed at the spatial infinity in this frame. Then the standard 
procedure is to perform the Fourier expansion of the initial data and follow the evolution 

^^The terms with gradients of the intrinsic curvature contribute with even higher derivatives of cj). 
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of each eigenmode separately. In this approach one indeed finds that in the boosted frame 
there are 2q eigenmodes. However, a straightforward analysis shows that only two of these 
modes have real frequencies. This means that the other modes grow exponentially either 
at positive or at negative times. In particular, if an observer living in the boosted frame 
is free to chose arbitrary initial conditions she will conclude that the system is unstable. 
However, the growing modes, while being legitimate solutions in the boosted frame, do not 
satisfy the condition ( 134|) in the original coordinate system. Thus in the latter system they 
are discarded as unphysical. 

The situation here is similar to the case of fields obeying second-order equation of motion 

- t;2A(^ = (35) 

with superluminal velocity, f > 1. The Cauchy problems for Eq. fl35l) are physically inequiv- 
alent when formulated in the original and highly boosted frames [T3| O US]. The latter 
corresponds to the case when the Cauchy slices intersect the future causal cone 

t = \x.\/v (36) 

of Eq. (1351) . On the other hand, all the Cauchy problems formulated on slices lying outside 
the cone fl36l) are equivalent. The difference between Eqs. fl33|) and fl35|) is that the former 
does not have a well-defined causal cone: the signal can propagate from the origin to any 
point at t > 0. Thus even arbitrarily small deviations from the original slicing qualitatively 
change the properties of the system. 

For what follows it is convenient to slightly generalize the above discussion. Consider the 
Lorentz covariant equation 

A^-'d.d,^ - (B'^'d^d^y ^ = 0, (37) 

where the symmetric matrices A^^", B^^ transform as tensors under Lorentz boosts. These 
matrices may depend on various fields present in the theory, in particular, they can depend 
on the field itself. Naively, Eq. (!33|) contains 2q time derivatives and thus describes q 
degrees of freedom. However this reasoning is not correct in general. The properties of the 
differential equation fl371) are characterized by the eigenvalues and eigenvectors of A^^ and 
B^'^ . The general study of the possible cases is beyond the scope of the present paper. Here 
we just point out the special case when B^^ has one hypersurface-orthogonal timelike eigen- 
vector with zero eigenvalue and 3 spacelike eigenvectors^ with non-zero eigenvalues. Then, 



^Note that because of symmetry of the matrix B^^'^ its eigenvectors are orthogonal. 
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in the frame defined by these eigenvectors the number of time derivatives in the operator 
B^^^d^dy is reduced. In general, in curved backgrounds or when the matrix B^^ is space-time 
dependent, the resulting operator still contains one time derivative, cf. Eq. fl38|) below. Still, 
the main conclusion is that the number of degrees of freedom described by an equation can 
be reduced compared to the naive expectation by proper choice of Cauchy slicing. 



Let us return to the Hofava gravity. Consider the higher derivative terms appearing in 
the equation of motion for the Stiickelberg field (f). To be concrete let us take the case of the 
action fl32l) . Then the term with four derivatives in the equation reads lf| 

This is precisely of the form fl57|) with B^^'^ = P^'^, q = 2. Thus we expect a reduction of the 
number of degrees of freedom in a certain frame. In the rest of this section we demonstrate 
that this is indeed the case for a large class of Hofava-type Lagrangians. Namely, we prove 
the following statement. Consider perturbations of the field around the background 0, 

= + X- (39) 

Then in the frame where the background is in unitary gauge, 

= t , (40) 

the linearized equation for x is first order in time derivative. 

For simplicity, we concentrate on the case when the Lagrangian for the Stiickelberg field 
depends only on first derivatives of the normal vector u^. Moreover, we assume that these 
derivatives enter into the Lagrangian through the extrinsic curvature /C^j^, 

£ = /:(m^/C^,). (41) 

This case covers all terms in the general Hofava-type Lagrangian except those involving 
spatial derivatives of the 3-dimensional curvature tensor Rijki- Indeed, the Gauss-Codazzi 
equation (!29|) implies that the terms polynomial in Rijki depend on only through the 
projector P^*^ and the extrinsic curvature. 



Explicitly in terms of 3+1 decomposition we have the operator identity 

P^^Vt.y,, = -A + ICu^Vx (38) 

valid for the action on scalar functions. 
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One observes that the equation of motion for the field has the form of the current 
conservation, 

V^J^ = , (42) 

where 



is the current related to the reparameterization symmetry ( |28l) . Let us demonstrate that the 
current is orthogonal to the gradient of (p, 

u^r = . (43) 

By explicit computation we find, 

This expression explicitly satisfies fl43p . Now we perform the separation of the field into the 
background and perturbations. The key observation is that in the frame defined by Eq. fHOl) 
the current contains exactly one time derivative of the linear perturbation x- Indeed, the 
perturbations of and the extrinsic curvature do not contain time derivatives of x- This 
follows from the explicit expressions 

Su" = -^Pf'^d^x , 



1 



/X 

where bar refers to the background values and 

is the proper acceleration of the congruency defined by u^. Then, by inspection of the 
expression (jUj) one finds that the only contribution with time derivative comes from the 
variation of the factor 1 / -\/X, 



To complete the argument we have to show that taking divergence of the current in the 
equation of motion ( H2l) does not bring more time derivatives. Due to the property ( l43l) one 
has 
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and thus (142!) can be written as 



Pl^s/^r - a,r = . 

When expanded to hnear order in perturbations, the first term contains only spatial deriva- 
tives of the perturbation of the current. Thus this term remains first derivative in time. 
Explicitly, the corresponding contribution reads, 

-v^v.x • 



Vx 

Another contribution with first time derivatives comes from the perturbation of a^y, 

Note that the two contributions are equal and sum up. This completes the proof. 

Two comments are in order. First, let us take a closer look at the equation of motion for 
the Stiickelberg field. For the sake of the argument let us consider the action (l32l) . Up to an 
irrelevant constant factor the current (jBl) reduces in this case to 

J'^ = ^P'^^V^/C . 

In the unitary gauge it takes the fornj^ 

jO = , / = N'^V'K . 

One observes that the equation of motion ( l42l) coincides with the secondary constraint ([9]). 
It is straightforward to check that the equation for is identical to the secondary constraint 
also for the general Hof ava type Lagrangian. This confirms the consistency of the Stiickelberg 
treatment. 

Second, the expression (HTj) is not the most general Lagrangian compatible with the 
reparameterization symmetry ( l28l) . Even if one restricts attention to Lagrangians with only 
first derivatives of u^, one can still add to fHTl) dependence on a'^ . In the unitary gauge the 
terms with translate into terms with spatial derivatives of the lapse, 

= , a} = V'N/N ; 
^^Recall that the covariant derivatives with Latin indices refer to the 3-dimensional metric jij. 
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such terms were not considered in As the hnear expansion of a*^ contains one time 
derivative of the perturbation x, the argument given above does not go through in this case, 
and the equation for the Stiickelberg field may contain more than one time derivativi^. In 
particular, adding the term (a'^)^ to the Lagrangian makes the equation second order in time. 

5 Instability and strong coupling 

We now perform a detailed study of the properties of the Stiickelberg field perturbations x- 
To make the analysis clear we restrict the study to the case when the Stiickelberg action 
contains only the term fl32|) . Moreover, we consider the decoupling limit, (1 — A) <^ 1. In 
this limit the backreaction of the Stiickelberg sector on the space-time geometry is negligible 
and one considers as propagating in a fixed background metric. At the same time this is 
precisely the limit where the theory is expected to approach GR. In this regime the non- 
linearities of the Stiickelberg field become important at energy scales much smaller than 
those of the other modes in the theory. This will allow to easily establish the strong coupling 
scale in this sector. 

We start by obtaining the quadratic action for the perturbations. According to the 
discussion of the previous section we work in the foliation defined by the background value 
of the field. Thus we fix the gauge (p = t, Ni = 0. Expanding the quantities in the action 
( 132|) to quadratic order one obtains 



This action explicitly reveals the properties of the extra mode discussed previously. Indeed, 
it contains only one time derivative of x, so that the resulting equation is first order in time. 
The term with time derivative vanishes whenever the gradient of the extrinsic curvature is 

This is readily understood in the unitary gauge, where the secondary constraint will now contain the 
time derivative of the lapse. 




(45) 



where we have used 



1C = K- 2V^NV'x - NAx + Nx^x + 2iVV*xViX 

- 

+ 2V,NxV\ + ATV.xV^X + —KViX'^'x - N^K^'VixV.x ■ 
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zero, and the field x becomes non-propagating. One also observes that due to the background 
equation of motion ^ this action is invariant under the shifts 



where ^(t) is an arbitrary function of time. This is recognized as the linearized form of the 
reparameterization (!28|) . This symmetry prevents x from having the ordinary x^ kinetic 
term. 

Let us analyze the equation of motion following from ( l45ll . As in Sec. [3] we are interested 
in the short wavelength limit, 

> 1 , 

where p is the momentum of the x-mode, and L is the characteristic length of the variation 
of the background. In this regime the dominant terms in the equation are those containing 
the largest number of derivatives of x- After using the equations for the background Q, the 
leading contributions to the equations of motion are 

2V'Kd,x + NA\ + 6V*iVV,Ax = , (46) 

where for future reference we retain the main subleading correction represented by the last 
term on the l.h.s. To find the solution of Eq. ( 146|) we use the same strategy as in Sec. [3l 
Namely, we restrict to the vicinity of a given point Xo- To the leading approximation the 
background fields in this vicinity can be considered as constant. One also assumes that the 
spatial metric 'jij at this point is fiat and the Cristoffel symbols vanish; this can be always 
achieved by performing a 3-dimensional diffeomorphism. Now it easy to obtain the leading 
behavior of the solution. One performs the Fourier expansion 

X oc e-'"*+'P" (47) 

and substitutes it into ( H6ll . Discarding the last - subleading - term one obtains the disper- 
sion relation 

^ = --^. (48) 

This coincides with the expressioiJ^ ( l23l) of Sec. [3] in the decoupling limit A ~ 1. 

We now discuss the first corrections in l/{pL) to the solution (H7!) . (HHI) . The reason for 
considering these corrections is that they qualitatively change the behaviour of the mode 



^^Recall that in the present section we work within the convention = 1. 
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making it either exponentially decaying or growing. As already pointed out in Sec. [3l in 
order to find the subleading corrections one has to go beyond the Fourier analysis and 
implement the WKB expansion in the vicinity of the point Xo- The details of this procedure 
are contained in the Appendix. Here we quote the main result. The eigenmode frequency 
acquires an imaginary part which is estimated as 

5uj ~ ip'^L . (49) 

The sign of this imaginary part depends on the direction of the mode momentum p relative 
to the gradients of the background. Those modes for which the imaginary part is positive are 
exponentially growing. Note that the rate of this growth is much faster than the characteristic 
background frequency 1/L. Thus we conclude that the Horava model suffers from fast 
instability at short scales. 

Another problem with the theory appears when one takes into account self-interaction of 
the field x- Let us first consider the case of fiat background. Then, computing the leading 
non-linear terms from the expansion of the covariant Lagrangian (152]) in the perturbation x 
we obtain, 

S = MM}^ J I (Ax)' + 2x ((Ax)' + 2d,xdiAx) } • (50) 

As expected, the quadratic part of the action does not contain time derivatives. Note, 
however, that the time derivatives do appear in the interaction. The form of this Lagrangian 
is restricted by the fact that x nonlinearly realizes the field-reparameterization symmetry 

(EHD 

The theory (150|) clearly has a dimensionful coupling constant (a/|1 — A| Mp)~^ which signals 
the presence of strong quantum coupling at high enough energies. The naive estimate of the 
strong coupling scale is provided by the inverse of this coupling [6], 

A„a»e = VI 1-A| Mp . (51) 

Note that this scale goes to zero in the putative GR limit A ^ 1. In particular, in this limit 
it is parametrically smaller than the "deep UV" scale Mp at which the higher-order terms 
of the Horava model could become important. 

We now argue that the scale of strong coupling for the action (l50!l is even lower than 
flSTj) . viz. zero. The physical reason is that due to the absence of time derivatives in the 
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quadratic part of the action, rapid fluctuations of the field x are not suppressed. Hence, the 
interaction terms with time derivatives blow up (see related discussion in [TB]). 

To make a quantitative statement let us regulate the action fl50|) by expanding on a nearly- 
flat but nontrivial background. This introduces the first order kinetic term for x as in fB^ . 
For momenta much larger than the scale defined by the background the leading order terms 
in the action are essentially the same as in flat space plus the background dependent kinetic 
term. Schematically, 

S = AL„e / {L-' v'xdiX + (Ax)' + xi^x)'} , 

where is the unit vector along the direction of the extrinsic curvature gradient, and L is 
the typical length scale of the background. The free part of the action is invariant under the 
scaling^ 

X I— s> b^^ X 
t^b'^t 
X^bx ■ 

Under this scaling the interaction term has dimension +4; thus it becomes relevant at short 
scales. To estimate the cutoff one performs the rescaling 

X ' ^ X ^ -^naiveX 5 

which brings the quadratic part of the action to the canonically normalized form. This yields, 

S= [d'x L'idd + (Ax)' + -^iiAxf 

.1 1 -^i-r 



'■naive 



where dot is now understood as the derivative with respect to the rescaled time i. From 
this expression one reads out the cutoff scales for the spatial momentum and the rescaled 
frequency: they are set by the appropriate powers of the unique coupling constant appearing 
in the interaction term. One obtains 

i\p u '■^naive ' 

A,. = L-'I'k^'L. . 



""^^Remarkably, we find here the same relative scahng of space and time as that proposed by Hofava for 
the deep UV to make the theory naively power-counting renormalizable. We have nothing more to say on 
that coincidence. 
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Using the relation Co = u/L'^ between the rescaled and the physical frequencies we obtain 
the physical frequency cutoff 

A. = L-'^'Alte ■ (52) 

Note that due to the non-relativistic structure of the theory the cutoff scales in spatial 
momentum and frequency are not equal. Instead, they satisfy the relation A^^ = A^L^, 
which is compatible with the dispersion relation of the extra mode. It is clear that both 
Ap and A^^ go to zero in the limit of flat background. Hence, the theory becomes strongly 
coupled at all scales. 

6 Hof ava gravity with projectable lapse as a ghost con- 
densate 

In this section we consider another version of Hof ava theory also proposed in [Ij. In this 
formulation the lapse function appearing in the action ([T]) is assumed to be projectable, 
which means that it does not depend on spatial coordinates, N = N{t). As we are about 
to see, the dynamics of the model in this case is very different from the case studied in the 
main body of the paper. We discuss it here only briefly, leaving a more thorough study for 
the future. 

In the projectable case the variation of the action with respect to N, instead of producing 
the local Hamiltonian constraint (jlj), gives the integral of Eq. (j4]) over the whole space. 
Such an integral constraint does not affect local physics. Thus, as far as local dynamics 
is concerned, the full set of equations of motion is providecjl^ by Eqs. ©-(ITj). Using the 
reparameterization of time one can set = 1 in these equations. Let us count the number 
of required Cauchy data. Out of 12 functions present in 'jij, 3 are constrained by Eq. ([5]). 
3 more functions are removed by residual spatial diffeomorphisms. Thus we are left with 
12 — 3 — 3 = 6 initial conditions. This is larger by 2 than in the GR case, suggesting the 
presence of an extra mode with second order evolution equation. 

The simplest way to study this new mode is by using the Stiickelberg formalism. First, 
notice that the constraint = 1 can be enforced by adding to the action ([1]) the term 

^P = /d^^v^Ar|(i,-l) , 

Clearly, Eq. ([8]) expressing the conservation of the Hamiltonian constraint is also absent. 



25 



where p is a Lagrange multiplier. From the relation (1301) . the generally covariant form of the 
previous action reads 

S, = J d'x^ P- (V^^V^^ - 1) . (53) 

Apart from this term the action for cj) contains higher-order contributions coming from the 
rest of the Horava Lagrangian. The latter contributions either contain more than two deriva- 
tives or describe non-minimal couplings to the metric, cf. Eq. (I3T11 . 

The theory fl53|) . together with higher order contributions, can be viewed as a special (non- 
minimally coupled) case of the ghost condensate model [TU]. The latter model is characterized 
by a Lagrangian 

Cgc = P(V ^(pV^cj)) + (terms with higher derivatives) , (54) 

where the function P{X) has a minimum at X = 1. This forces the field to develop a 
non-zero time- like gradient. The Lagrangian (153|1 is the "sigma model limit" of (!54|) when 
the function P gets replaced by the constraint. Note, however, an important difference in 
the interpretation of the higher derivative terms in the ghost condensate model and Hofava 
gravity. In the standard approach to the ghost condensate they are considered as higher 
order terms in the effective field theory expansion, while according to Horava they should 
be taken at face value and determine the UV properties of the system. 

Let us get more insight into the dynamics of the field described by (!53l) . For the sake 
of the argument we omit the higher order terms. This would correspond to the IR limit of 
the putative UV complete theory. However, we will see shortly that the 0-sector exhibits a 
number of pathologies which make the UV completion problematic. 

The action yields the following equations of motion, 

V,(PW(P = 1 , (55) 
V^(pV'^0) = . 

These equations are equivalent to the equations of motion of an ideal pressureless fluid (dust) 
with energy-momentum tensor 

Tf,u = P Vf,(f)Vu(j) . 

We see that p and are identified respectively as the energy density and velocity of the 
effective fluid. This result is in agreement with the findings of [TT] where it was proposed to 
interpret the effective dust-like component arising in the projectable version of the Hofava 
gravity as dark matter. However, we point out that a well-known property of pressureless 



26 



fluids is that they develop caustics, i.e. there are space-time regions where the fluid velocity 
is ill-defined (see p[7] for discussion of this topic in the context of field theory models of dark 
matter). The formation of caustics is easy to understand in the decoupling limit when the 
backreaction of the fluid on the geometry is negligible. Then the fluid particles move along 
geodesies without feeling each other. Given a general inhomogeneous initial distribution of 
particle velocities, their trajectories will cross forming a caustic. While this process does 
not pose a problem for real dust where it leads to virialization, the formation of caustics 
means an inconsistency in the case of a scalar theory, as the field is not differentiable at the 
caustics. Thus, we conclude that for generic initial configurations the theory described by 
the action flS5]) breaks down after finite time evolution. Note that the formation of caustics is 
a general problem of the ghost condensate action (15^ |18j . In that context, it was suggested 
[T8] that the problem might be resolved by the effect of higher derivative term which make 
the fluid particles deviate from the geodesic motion. This hope is absent in the special case 
of the action (l53ll : irrespective of the higher-order terms the fluid particles move exactly 
along geodesies as long as the constraint is present. 

The consistency of the Lagrangian flS^ is also challenged at the quantum level. To get 
the flavor of the problem let us again omit the higher order terms and proceed to quantize 
the model in the canonical formalism. For simplicity we also neglect the backreaction of the 
field on the metric and assume the latter to be flat. The canonically conjugate momenta 
for the variables 0, p are 

TTp = . (56) 
The constraint ( l55l) in canonical variables takes the form 

4 - (9^<l>r = 1 . (57) 
P 

Equations ( l56i) and ( l55i) form a pair of second class constraints and enable to eliminate the 
variables p, it p. In this way one obtains the Hamiltonian, 

H = 7r^Vl + (9,0)2 . (58) 

Note that in our case the elimination of constraints does not modify the Poisson brackets of 
the remaining variables, 

{7r^(x), 0(y)} = 5(x - y) , etc. 
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The quantization of the theory proceeds now in the standard way by imposing the canonical 
commutation relations on vr^ and 0. 

The Hamiltonian (|58|) certainly looks unusual. To get insight into its properties, let us 
expand it around the background = t, vr^ = po- Thus we write 

(t^ = t + xl a/Po , vr,^ = Po + T^x\fPo 

and obtain 

n = \{daf + ^^Adaf + ... . (59) 

As an example we included here one of the interaction terms. Now, ( l59l) describes a theory 
with dimensionful coupling This implies that the theory gets strongly coupled at 

the cutoff scale A < p^^ . For the dark matter interpretation of the 0-sector, taking the 
present-day value for the average density of dark matter, one obtains the cutoff 

A < 10-=^eV , 

which is unacceptably low for a candidate theory of quantum gravity. In fact, the cutoff 
in the theory fl59|) may be even lower, viz. zero. This can be argued from the fact that 
the quadratic part of the Hamiltonian (!59|) does not contain the momentum vr^. Thus the 
quantum fluctuations of tt^ are not suppressed and the interaction terms containing tt^ blow 
up. The careful analysis of this issue is beyond the scope of this papec^. 

To summarize, we found that in the decoupling limit, when the backreaction of the extra 
mode on the metric can be neglected, the projectable version of the Hofava gravity suffers 
from caustic and low cutoff problems. This suggests the inconsistency of the model in its 
present form. 



7 Conclusions and discussion 

In this paper we have studied Horava's proposal for quantum gravity pp from the low-energy 
perspective. We have mainly concentrated on the "non-projectable" version of the model. 
We have uncovered the additional scalar degree of freedom arising from the explicit breaking 

""^^Note that the canonical transformation from the variables (j), tt^ to t^x time-dependent. Taking this 
properly into account eliminates the term linear in vr^ in the Hamiltonian. 

Another issue which we do not address in this paper is the effect of the higher order terms on the 
power-counting. 
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of the general covariance and analyzed its properties in detail. A peculiarity of the new mode 
is that it satisfies an equation of motion that is of first order in time derivatives; this means 
that it adds just one direction to the phase space, or 'half a degree of freedom. At linear level 
the mode is manifest only on spatially inhomogeneous and time-dependent backgrounds. We 
have demonstrated the existence of the extra mode in two physically equivalent ways. First 
we performed the analysis in the original ADM-like formulation of Horava and identified 
the mode among the perturbations of the metric. Second, we made use of the Stiickelberg 
formalism which amounted to restoring the general covariance of the model. This required 
introducing a scalar field describing the foliation structure of the Horava model. It is worth 
stressing that this procedure did not add any new physical degrees of freedom compared 
to the original formulation because the scalar field was introduced simultaneously with a 
new local symmetry. The latter gauge symmetry is the part of general covariance which is 
explicitly broken in the Hofava's original formulation. The Stiickelberg approach allowed us 
to transfer the extra mode into the fluctuations of the foliation structure and to perform a 
detailed study of its properties. 

We found two serious problems associated with this mode. First, the mode develops 
very fast exponential instabilities at short distances. Second, it becomes strongly coupled at 
extremely low cutoff scale. Due to the non-relativistic nature of the theory the cutoff scales 
in spatial momentum and frequency are different. They both depend on the curvature of 
the background metric and go to zero when this curvature vanishes. These features allow to 
conclude that Hofava's proposal is inconsistent in the present form. 

We have also discussed the "projectable" version of the Horava model. We have argued 
that this version can be understood as a certain limit of ghost condensation. In this case, 
the theory propagates a whole degree of freedom with second order equation in time for any 
background. In this sense the projectable version is better behaved than the non-projectable 
one. However, the preliminary analysis of the model in the regime when the backreaction 
of the extra mode on the space-time geometry can be neglected suggests that the model is 
still problematic since the additional fleld generically forms caustics and, again, has very low 
strong coupling scale. 

Let us comment on the possible ways in which the problems of the theory can be ad- 
dressed. The comparison of projectable and non-projectable cases suggests that the only way 
to make the extra mode well-behaved is to promote it to a full-fledged scalar. There are two 
strategies to do so. One possibility is to relax the requirement of invariance under reparam- 
eterizations of time and reduce the symmetry group of the theory down to (time-dependent) 
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spatial diffeomorphisms. As a result one obtains some version of the ghost condensate model 
[TU] which is known to be a consistent effective theory up to relatively high scales. The 
problem with caustics can be also addressed in this framework [18]. In the unitary gauge 
the theory would be described by a generalization of Hofava action containing arbitrary 
functions of the lapse N. Needless to say, it is unclear if the appealing UV properties of the 
Hofava proposal can be preserved within this approach. An argument suggesting that it may 
be impossible is the break-down of the black hole thermodynamics in the ghost condensate 
model [T9l [20] . and more generally, in theories with low-energy violation of Lorentz invari- 
ance [21]. Indeed, the break-down of thermodynamics indicates violation of unitarity in the 
underlying theory, which makes construction of a renormalizable and hence UV complete 
quantum theory having the ghost condensate as its low-energy limit problematic. 

Alternatively, one can promote the extra mode to a whole scalar without breaking the 
invariance under foliation-preserving diffeomorphisms. To achieve this one has to add to 
Hofava's Lagrangian terms that are non-linear in and respect the symmetry. An example 
of such term is A^~^(VjA^)^. In the Stiickelberg picture around fiat background this corre- 
sponds to the addition of the term (dix)^. In this case also the equation for the Stiickelberg 
field X becomes of second order in time for any background. Still, the terms added to the 
Lagrangian in the Stiickelberg language are higher order in total number of derivatives. It 
remains an open question if this feature does not lead to any pathologies. The problems 
with the break-down of black hole thermodynamics would arise in this approach as well: the 
field X generically has different propagation velocity compared to the graviton (helicity-2 
excitation), thus violating the Lorentz invariance at low energies. It is unclear at present if 
and how these problems can be resolved. 

Finally, we mention another interesting outcome of our study. We observed that the 
equation of motion for the Stiickelberg field in Hofava gravity has a peculiar structure. 
This equation is explicitly higher order in covariant derivatives. Consequently in a general 
coordinate frame it is higher order in time. Generically, this would imply the presence of 
additional ghost modes. However, this reasoning is incorrect in our case: there is a unique 
preferred frame where the number of time derivatives in the equation is drastically reduced. 
Solving the Cauchy problem in this frame requires less initial data than in a general frame. 
The difference can be traced to the intrinsically non-relativistic nature of the theory which 
implies that the two frames are physically inequivalent. Technically this manifests itself in the 
fact that the boundary conditions at spatial infinity are not equivalent in the two frames. 
We pointed out that the reduction of degrees of freedom in a preferred Cauchy slicing is 
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generic for a wide class of higher order covariant derivative operators. This opens up the 
possibihty to construct a new class of consistent higher derivative theories with equations of 
motion based on these operators. We leave the investigation of these issues for the future. 
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A WKB expansion for the extra mode 

In this Appendix we construct the solution of Eq. f H6l) using the WKB method. The ap- 
plication of this method is legitimate in the case when the frequency uj and momentum p 
of the field x ^i-re much larger than where L is the characteristic space-time scale of 
the background. The leading-order approximation to the solution is obtained in the main 
text and is given by the plain wave P7|) with the frequency related to the momentum by 
the dispersion relation ( HHl) . Our aim here is to obtain the order l/{pL) corrections to this 
solution. 

We work locally in the vicinity of the point Xq, which without loss of generality we assume 
to coincide with the origin of the coordinate frame. By the 3-dimensional diffeomorphism we 
can make the spatial metric at this point flat, and its first derivatives vanish. This implies 
that the spatial covariant derivatives in Eq. (l46l) can be replaced by ordinary ones in our 
approximation. We now use the following ansatz for the field x- 



X oc exp 



- iuj^^h + ipix' - idojt + i^t^ + ibitx" + ^Cijx'x^ . (60) 



Here u;^^^ is the leading-order frequency fHSl) evaluated at Xo] Su is the frequency shift; and 
(2 , 62 5 c^ij are the coefficients in the Taylor expansion of the WKB phase at Xo- Note that if 
we want the ansatz fl60l) to represent a small correction to the plain wave solution fH7|) in 
the region t, -C L, we must require 

Suj < , a ~ 00^^^ L-^ , bi, dj ~ pL-^ . (61) 
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One proceeds by evaluating the derivatives of x appearing in Eq. (146|) . Keeping the first 
subleading corrections one obtains, 

diX = + Pi^^ - PM - PibjX^ + uj'-^^bit + ijj^^'^CijX^ + ^h^] x , (62) 

^^X = [p^ + ^P^PiH + Ap^PiCijX^ - AipiPjdj - 2ip^Cii] x ■ (63) 

The last term in Eq. (H6l) is already subleading because of the additional derivative of the 
background. Thus, it is enough to evaluate the corresponding x-derivative to the leading 
order 

diAx = -ipiP^X ■ (64) 

Finally, the inhomogeneity of the background is taken into account by expanding the coeffi- 
cients in (1461) in a Taylor series, 

diK = diKo + didjKoX^ + diRot , (65) 
N = No + diNox' + Not , (66) 

where the quantities with the subscript "o" are evaluated at the origin Xo- We will omit this 
index in what follows. 

The next step is to substitute the expressions (162|) - (166|) into Eq. (H6|l and require the 
vanishing of the Taylor coefficients up to linear order in coordinates. This yields the following 



system of equations, 

- 2pidiKa + (2cu(°)9ii? + ANp'^pi)bi + 2j°^pidiR + p'^N = , (67) 

- 2pidiKbj + {2uj^^^diK + ANp^pi) aj + 2j^^p,didjK + p^djN = , (68) 
pAKSu; + idiKbi - iN{2piPjC,j + ph,i) - Sip^pAN = . (69) 



Note that within the assumptions (I6T1) the first term in Eq. (l68ll is much smaller than the 
rest and we can neglect it. From this equation, the solution for Cij reads: 

where 

Ui = 2uj^^AiK + Ap^piN , 
Vj = 2J^^p'didjK + p'^djN . 



32 



It is easy to check that Cij satisfies the estimate (!6T|) . We do not need to solve the system 
in the full generality: any special solution is sufficient for our purposes. Thus we 
set bi = and find a from fl67|) . Finally, from Eq. fl69|) we obtain the frequency shift 



This expression is purely imaginary. Analyzing the order of magnitude of the terms entering 
into it one obtains the estimate (H9ll used in the main text. 
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